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Abstract 

We study the radiative pion decay tt + — > e + v e ^ within nonlocal chiral quark models that include 
wave function renormalization. In this framework we analyze the momentum dependence of the 
vector form factor Fv(q 2 ), and the slope of the axial- vector form factor FA(q 2 ) at threshold. Our 
results are compared with available experimental information and with the predictions given by 
the NJL model. In addition we calculate the low energy constants and £q, comparing our results 
with the values obtained in chiral perturbation theory. 

PACS numbers: 12.39.Ki, 11.30.Rd, 13.20.Cz 
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I. INTRODUCTION 



The radiative pion decay 7r + — > e + v e r ) is a very interesting process from different points 
of view. According to the standard description, the corresponding decay amplitude consists 
of the inner bremsstrahlung (IB) and structure-dependent (SD) terms. The former can be 
associated with the diagrams in which the photon is radiated by the electrically charged 
external legs (either pion or lepton), while the SD terms correspond to the photon emission 
from intermediate states generated by strong interactions. Since the IB contribution turns 
out to be helicity suppressed, this process is an adequate channel to study the SD amplitude, 
which provides information about strong interactions in the nonperturbative regime. 

The SD contribution can be parameterized through the introduction of vector and axial- 
vector form factors, Fy(q 2 ) and Fa{q 2 ) respectively, where q 2 is the squared invariant mass 
of the e + z/ e pair From the experimental point of view, recent analyses |2J of n + — » e + z/ e 7 
decays have allowed the measurement of -FV(O), -Fa(O), and the slope of the form factor 
F v (q 2 ) at q 2 = 0. In addition, ongoing experiments are expected to reach enough statistics 
to determine the slope of F^iq 2 ) at q 2 = in the near future. On the theoretical side, the 



analysis of the 



brm factors has been carried out in the framework of Chiral Perturbation 



Theory (%PT) 3] and effective meson Lagrangian methods [4], which provide a good de- 
scription of low-energy meson phenomenology. However, one can also address the question 
of how the form factors are connected to the underlying quark structure. Due to the nonper- 
turbative nature of the quark-gluon interactions in the low-energy domain, to address this 
issue one is forced to deal with models that treat quark interactions in some effective way. 
In this sense, the widely studied Nambu-Jona-Lasinio (NJL) model is the most popular 
schematic quark effective theory for QCD 5^7]. In this model quarks interact through a 
local, chiral invariant four-fermion coupling. Now, as a step towards a more adequate de- 
scription of QCD, some extensions of the NJL that include nonlocal interactions have been 
proposed in the literature (see Ref. js] and references therein). In fact, nonlocality arises as 
a natural feature of several well established approaches to low energy quark dynamics, such 
as the instanton liquid model J9jl, the Schwinger-Dyson resummation techniques 10(, and 
also lattice QCD calculations [ljj-tl3|. In the last years nonlocal chiral quark models have 
been applied to study different hadron observables with significant success 14j-|22|. 

In a previous work (23J we have addressed the analysis of the vector and axial-vector 
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form factors in 7r + — > e + u e 'j in the framework of nonlocal chiral quark models. We have 
determined the values of i*V(0), -^4(0) and the slope of Fy(q 2 ) at q 2 = for different 
parameterization sets, comparing the results with the experimental measurements and the 
corresponding values obtained in the NJL model. In this paper we extend these previous 
results, analyzing the vector form factor for values of q 2 up to 1 GeV 2 , and studying the 
slope of FA{q 2 ) at threshold. While the latter is interesting in view of the comparison 
with future experimental results, the study of the behavior of Fy(q 2 ) presents relevant 
theoretical motivations. In fact, in the isospin limit, and assuming the conserved- vector- 
current (CVC) hypothesis, the vector form factor can be directly related to the form factor 
F™ 71 * (q 2 ) associated with the vertex ^77*, where q 2 is the photon virtuality in Euclidean 
space. Thus our results can be compared with the experimental information on F 7171 * (q 2 ) , 
taken e.g. from processes as n~p — > 7T°n or e + e~ — > e + e~7r°, and the compatibility with the 
usually assumed single-pole behavior can be analyzed. The CVC hypothesis is supported by 
the measured value of -FV(O), which turns out to be fully consistent with the corresponding 
value taken from F 71 " 77 * (q 2 ) at q 2 = [24]. We notice that the calculation of Fy(q 2 ) in the 
framework of nonlocal models turns out to be technically involved, due to the presence of 
poles and cuts arising from the dressed quark propagators. Finally, to complete our analysis 
we also study the compatibility of our nonlocal models with effective meson theories such 
as xPT. In fact, it is seen that Fa(0) and the pion charge radius (r 2 ) n are related to the 
so-called %PT low energy constants (LECs) £5 and £q, which encode the information of 
quark and gluon degrees of freedom. In %PT these LECs are input parameters taken from 
phenomenology, thus the compatibility with the predictions arising from effective quark 
models is a measure of the ability of these models to account for the underlying strong 
interaction physics. 

The article is organized as follows: in Sect. II we describe the models and quote the 
analytical results for the vector form factor Fv{q 2 )- The numerical results for vector and 
axial-vector form factors are presented and discussed in Sect. III. In Sect. IV we quote the 
values for the LECs obtained within our models, while in Sect. V we present our conclusions. 
The full analytical expressions for the axial- vector form factor are quoted in Appendix I. 
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II. FORMALISM 



The amplitude for the process tt + — > e + u e (q) + 'y(k) in Minkowski space can be written 
as Q 



G 

M = —= e cos 8 C 
V2 



V2 f v \(q + k) a L ait - l v 



(q-k) 



+l v \-ie^ k a qt 



+ [%K - (q ■ k)} 



F A (q 2 



(1) 



with (q + k) 2 = m 2 , k 2 = 0. Here Gp and 9c stand for the Fermi constant and the 
Cabibbo angle, respectively; e M is the photon polarization vector, Z M is the lepton current, 
L a/1 is a lepton tensor, and Fy(q 2 ) and FA{q 2 ) denote the vector and axial- vector hadronic 
form factors mentioned in the Introduction. In this work we are interested in the study of 
these form factors in the context of nonlocal SU(2) chiral models that include wave function 
renormalization. These models are defined by the following Euclidean action 2l|, [22] : 



St 



d 4 x <ijj(x)(—i^ + m c ) , ijj(x) 



Gs 
2 



ja(x)j a (x) +jp(x)j P (x) 



(2) 



Here m c is the current quark mass, which is assumed to be equal for u and d quarks, while 
the nonlocal currents j a (x), jp(x) are given by 



Ja[X) 



d z g{z) ^[x + - ) T a ip(x- - 



jp( x ) = / d4z f( z ) i> ( x + 



ip [x — 



(3) 



where T a = (l,i7 5 r) and u(x') d v(x) = u(x')d x v(x) — d x >u(x')v(x). The functions g(z) and 
f(z) in Eq. ([3]) are nonlocal covariant vertex form factors characterizing the corresponding 
interactions. In what follows it is convenient to Fourier transform g(z) and f(z) into momen- 
tum space. Note that Lorentz invariance implies that the corresponding Fourier transforms, 
g p and f p , can only be functions of p 2 . 

In order to deal with meson degrees of freedom, one can perform a standard bosoniza- 
tion of the theory. This is done by considering the corresponding partition function 
Z = f UipVip exp[— Se], and introducing auxiliary fields cr 2 (x), n(x), where a\^{x) 

and 7r(x) are scalar and pseudoscalar mesons, respectively. An effective bosonized action is 
obtained once the fermion fields are integrated out. To treat that bosonic action we assume, 
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as customary, that a± t 2 fields have nontrivial translational invariant mean field values, while 
the mean field values of pseudoscalar fields 7Tj are zero. Thus we write 



ai(x) — a\ + 5ai(x) 



a 2 (x) = k p a 2 + 5a 2 (x) 



-k(x) = Snix) 



(4) 



Replacing in the bosonized effective action and expanding in powers of meson fluctuations 
we get 



S 



bos 



s 



MFA 



+ sr d + 



Here the mean field action per unit volume reads 



cMFA 



2G 



' {al + >4ai)-AN c [ £^ h, 



(27T) 



with 



(5) 



(6) 



(7) 



z p = (1 - cr 2 f p ) 1 , m p = z p (m c + cti g p 

The minimization of S|? FA with respect to leads to the corresponding "gap equations". 
The quadratic terms in Eq. (|3j) can be written as 

1 r d 4 p 



oquad 



J2 G M (p 2 ) 5M(p) 5M(-p) 



(8) 



where a and a' fields are scalar meson mass eigenstates, defined in such a way that there is 



no o — o"' mixing at the level of the quadratic action. The explicit expressions o 



G M (p 2 ), 

as well as those of the gap equations mentioned above, can be found in Ref. 22]. Meson 
masses can be obtained by solving the equations Gm{— t^m) — 0, while on-shell meson-quark 
coupling constants GMqq are given by 

dG M (p 2 



(G 



Mqg) 



dp 2 



(9) 



o2 



As in Ref. 



22], we will consider here different functional dependencies for the vertex form 



factors g p and f p . First, we consider a relatively simple case in which there is no wave 
function renormalization of the quark propagator, i.e. f p — 0, z p — 1, and we take an often 
used exponential parameterization for g p , 



g p = exp (-p 2 / A o) • 



(10) 



The model parameters m c , Gs and A are determined by fitting the pion mass and decay 
constant to their empirical values m n = 139 MeV and /„- = 92.4 MeV, and fixing the 



chiral condensate to the phenomenologically acceptable value (qq) 1 ^ 3 = —240 MeV. In what 
follows we refer to this choice of model parameters as set A. Second, we consider a more 
general case that includes the jp{x)jp(x) interaction, i.e. we consider a nonzero wave function 
renormalization of the quark propagator. We keep the exponential shape in Eq. ffTUj) for the 
form factor g p and assume also an exponential form for f p , namely 

f p = exp (V/A?) • (11) 

Note that the range (in momentum space) of the nonlocality in each channel is determined 
by the parameters Ao and Ai, respectively. The model includes now two additional free 
parameters, namely Ai and k p . We determine here the five model parameters so as to 
obtain the desired values of m n , f w and (qq) 1 ^ 3 , and in addition we fix m c = 5.7 MeV and 
impose the condition zJO) = 0.7, which is within the range of values suggested by recent 



lattice calculations 
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131 ] . We will refer to this choice of model parameters and form 
factors as parameterization set B. Finally, we consider a different functional form for the 
form factors, given by 

= 1 + a z am fmjp) - m a z f z (p) _ 1 + a z 

9p l + a z f z (p) a m -m a z p 1 + a z f z (p) ' v ' 



where 



Up) = i + (p 2 /a 2 o) 3/2 ■ Up) = [1 + ^/A?)]- tt/a . (13) 



-5/2 



With this parameterization, the wave function renormalization of the quark propagator 
and the effective quark mass have the simple expressions Z p = 1 + a z f z (p) and m p = 
m c + oc m f m {p). As shown in Ref. |22j, taking m c = 2.37 MeV, a m = 309 MeV, a z = —0.3, 
A = 850 MeV and Aj = 1400 MeV one can very well reproduce the momentum dependence 
of mass and renormalization functions obtained in lattice calculations, as well as the physical 
values of m n and /„-. We will refer to this choice of model parameters as parameterization 
set C. The parameter values for all three parameter sets, as well as the corresponding 
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predictions for several meson properties, can be found in Ref. 

In order to derive the form factors we are interested in, one should "gauge" the effective 
action Se by introducing the electromagnetic field A^(x) and the charged weak fields WMx). 
For a local theory this "gauging" procedure is usually done by performing the replacement 

d^d^+iG^x) , (14) 
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where 



(15) 



2 V2 

with gj^/(8M^) = Gp/v / 2andr ± = (r 1 ±zr 2 )/2. In the present case — owing to the nonlo- 
cality of the involved fields — one has to perform additional replacements in the interaction 
terms, namely 



ip(x-z/2) -> W G {x,x - z/2) ip(x-z/2) 
ip\x + z/2) -> ^(x + z/2) W G (x + z/2,a;) 



(16) 



Here x and z are the variables appearing in the definitions of the nonlocal currents [see 
Eq.flH])], and the function Wa(x,y) is defined by 



W G (x,y) = P exp 



(17) 



where r runs over an arbitrary path connecting x with y. Once the gauged effective action 
is built, the explicit expressions for the vector and axial- vector form factors can be obtained 
by expanding to leading order in the product 5tt + Wj~. 

In the case of the vector form factor, one gets only a contribution associated with the 
triangle diagram represented in Fig. la. This is given by 



Fv{q 
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where for convenience we have used the shorthand notation = £ 2 + m 2 and if = £±r/2. 
The terms t\ and ti in Eq. ffl8|) read 
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The integrals in a h (£, q) have their origin in the path dependence introduced through the 
"gauged" nonlocal effective action. The result in Eqs. (I18ti20p corresponds to the use of a 
straight line between points x and y in Eq. (IT7|) . 

For the case of the axial form factor the calculation is more involved since it receives 
not only a contribution from the triangle diagram in Fig. la (as occurs in the local NJL 
model) but also from other diagrams, which are represented in Figs, lb-le. Since the result- 
ing expressions are rather long, we choose to relegate them to Appendix I where they are 
presented in detail. We recall that we have been working in Euclidean space, therefore our 
expressions for both Fy(q 2 ) and F^q 2 ) correspond to Euclidean q 2 . In order to obtain the 
form factors defined through the tt + — > e + u e ^ amplitude in Minkowski space [c.f. Eq. 
one should change q 2 — > —q 2 . 

III. NUMERICAL RESULTS FOR THE VECTOR AND AXIAL FORM FACTORS 

In this section we present and discuss our numerical results for the vector and axial- 
vector form factors. Let us start with the vector form factor Fv{q 2 )- We will only restrict 
to values of q 2 which are in the region where one can rely on the applicability of effective 
chiral models, i.e. q 2 < 1 GeV 2 . It is important to stress that, even within that region, the 
numerical evaluation of Eq. ffTSj) presents some technical difficulties due to the presence of 
poles in the corresponding integrand. Namely, once q 2 increases beyond a certain threshold 
(about 0.1 — 0.2 GeV 2 for the parameterizations considered here) the functions Dg + k and 
Di- q in Eq. (1181) can become zero for real values of £4. Thus, in order to properly evaluate 
the integrals, the corresponding residues have to be conveniently added. Since for non-local 
interactions the quark propagators might have several different poles, this procedure has 
to be repeated each time a new pole threshold is overpassed. In the case of set C the 
situation is even more complicated due to the existence of a cut arising from the particular 
q 2 dependence of the interaction form factors, Eqs. (fT2l) and f|T3|) . In any case, once these 
issues are taken into account, a smooth q 2 dependence of the vector form factors is obtained. 

Before discussing in detail our predictions for the form factors associated with the charged 
pion radiative weak decay, let us note that one can also consider the related decay processes 
7T° — > 77 and 7T° e + e~7. For these processes the amplitude involves the 7r°77 vertex 
form factor F 71 " 77 * (q 2 ), where q 2 stands for the squared invariant mass of the virtual photon, 
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defined here in Euclidean space. Experimental results can be well described by a single-pole 
parameterization, 

77^77* (r\\ 

Here it is important to mention that, in the chiral limit, the U(l) axial anomaly leads to 
the constraint F n ^' (0) = l/(47r 2 / 7 r r L^ 0.274 GeV~\ which is well satisfied by the non-local 
models of the type considered here [161 ] . Moreover, assuming the validity of the CVC hypoth- 
esis together with isospin arguments it is possible to prove that F v (0) = m T+ F n ^*(0)/V2. 
In the case of the non-local models under consideration the validity of this relation also 
extends to finite values of q 2 , namely 

F v (q 2 )= m7T+ F^*(q 2 )/V2. (22) 

The reason for this is that, in addition to respect CVC and isospin symmetry, the models do 
not include explicit vector channels in the quark-antiquark current- current interactions. As 
discussed below, differences between the isoscalar and isovector vector channel interactions 
could lead to violations of the relation in Eq. ( f22|) at finite q 2 . 

Our results for Fy(q 2 ) are displayed in Fig. 2, where we plot the curves corresponding 
to the form factors introduced in the previous section. In addition, in Fig. 2 we show the 
results obtained from the combination of Eq. (1221) with the parameterization in Eq. (|2~T|) . The 
parameters are taken from the experimental fit to 7r° — > 77* data, which yields F 71 " 77 * (0) = 
0.284(8) GeV -1 , a' v = 0.032(4) [24|. In general, it is seen that the results for our parameter 
sets A, B and C turn out to be similar in the whole range of q 2 values considered. The 
particular values at q 2 = are given in Table I. As already remarked in Ref. [23], one expect 
these values to be basically coincident given the above mentioned constraint in the chiral 
limit. Concerning the q 2 dependence, it is interesting to note that for all three parameter 
sets the behavior can be very well reproduced by a single-pole fit of the form 

M?) = . / v( °' -, ■ (23) 
1 + a v q z /m* 

In view of the relation Eq. (1221) . in our model ay can be strictly identified with the parameter 
a' v introduced in Eq. (12Tj) . 

On the other hand, it is possible to analyze the q 2 dependence of the vector form factor 
considering the experimental data obtained from the process 7r + — > e + v e ^. For low values 
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of q 2 one can carry out an expansion of the type 



F v {q 2 ) = F V (Q) 



q 2 f a 4 
1 — Ay — 2 1" -V — T 

m+ m 4 



(24) 



where the minus sign arises from the definition of Ay = m 2 Fy(0)/Fy(0) for q 2 in Minkowski 
space. In Table I the values of Ay, Ay and ay are given. We observe that ay is slightly 
different from the value of Ay. The discrepancy can be taken as measure how much the 
actual q 2 dependence of the vector form factors deviates from the single-pole behavior. This 
slight deviation can be also seen from the closeness between Ay and Ay, which should be 
equal in the single-pole limit. In general, it is seen that for all three parameterizations the 
results for Ay and ay are lower than experimental estimates (notice that some enhancement 
is found when going from sets A, B to the lattice-inspired set C). We remark that, given the 
fact that the calculated values of Fy(q 2 ) are well adjusted by the single-pole fit in Eq. ( f23i) . 
the discrepancies between the model predictions and the empirical fit observed in Fig. 2 
can be traced back to the values of ay (or Ay), i.e. the slopes of the curves at q 2 = 0. 
In fact, these discrepancies are not unexpected, since our model does not include vector- 
vector interactions. The magnitude of the corresponding corrections can be estimated by 



considering e.g. the extended NJL model studied in Ref. [25]: indeed, it is seen that the 
vector contribution to 7r° — > 77* has the right order of magnitude and sign in order to 
account for the discrepancies. It is also interesting to point out that the contribution from 
the vector channel could be different for ay and a' v , even when isospin symmetry is preserved. 
This would be achieved if one has different interactions in the vector-isoscalar channel and 
in the vector-isovector channel. As a final comment concerning the vector form factor, let us 
discuss the path dependence of our results. Contrary to what happens at q 2 = 0, the values 
obtained at finite q 2 depend on the path considered in Eq. ( TT71) . As usual, we have chosen 
a straight line path for the calculations of the path dependent quantities a^(£, q), Eq. fT20l) . 
To estimate the significance of this path dependence we have performed single-pole fits to 
the predictions obtained by neglecting the contribution of the the corresponding terms in 
Eq. ( 1191) . The resulting values for ay turn out to differ from the values listed in Table I by 
less than 0.5 %. 

We turn now to the axial form factor F^g 2 ). Contrary to the case of the vector form 
factor, FA^q 2 ) is not experimentally accessible for (Euclidean) q 2 > 0. As stated, the value 
of Fa{0) has been measured from n + — > e + v e ^ decay, and forthcoming experiments offer 
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the possibility of determining the slope -^(0) = —Fa(0) Aa/^+- The model predictions 



for these quantities are given in Table II. As it was discussed in Ref. [23j, for Fa(0) the 
triangle diagram in Fig. la turns out to be the dominant one, giving at least 98% of the 
total value for all three parameter sets considered here. In the case of the slope the triangle 
diagram is still the dominant one, but the contributions from diagrams in Figs, lb and le 
are also significant. One should note that nonlocal models of the type considered here lead 
to values of Fa(0) which are significantly different from those of -FV(O). This is remarkable, 
since other approaches like the NJL model 26( or the spectral quark model 27j tend to 
give -FV(O) ~ Fa(0). Given the triangle diagram dominance mentioned above, the origin 
of this difference can be traced back to the different dressing of the 7 M and 7^75 terms in 



the coupling of the W to the quarks 



22j. From Table II it is seen that the predictions for 



Fa(0) in nonlocal models are significantly closer than the experimental value than the result 
obtained in the standard NJL model. Regarding the slope parameter A^, we find that our 
predictions are quite similar for the three parameter sets considered, and turn out to be 
about 3/5 smaller than the value obtained in a resonance effective model 4\. 



IV. PREDICTIONS FOR LECS OF CHIRAL PERTURBATION THEORY 



In the previous section we have focused our attention on the ability of our quark models 
to reproduce the main features of the vector and axial-vector form factors. An alternative 
point of view (see for example Ref. j^sj) is to consider quark models as the generators of the 



pion xPT Lagrangian 



29]. ;\RT describes the low energy physics of pions in a universal way, 



once the order in the momentum and chiral breaking expansion (i.e. the order in the chiral 
expansion) is specified. Different scenarios for quark models will lead to Lagrangians 
with different values of the corresponding low energy constants (LECs). In this section we 
analyze the connection between our quark scenarios and the ;^PT Lagrangian up to the 
fourth order in the chiral expansion. 

The pionic Lagrangian that follows from performing a fourth order gradient expansion 
of our model bosonized action in the presence of external fields can be obtained using, for 
example, the method described in Ref. 28J]. On general grounds, it should read 



jC — £-2 + £4 , 



(25) 
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where 



c2 



£2 = y V M f/ T V»U + f X T U , (26) 



£4 = h (v M t/ T v^)' + £ 2 {v^u T v u u) {v»u T v v u) + 4 { x T uy + 

4 (V M x T V^f/) 2 + 4 U T F^F^U + 4 WU T F^VU + ... , (27) 



Here, {7 = (s,pi) is a four-component real 0(4) vector field of unit length, [/ T [7 = 1, and 
X = 2B(S,Pi) is another 0(4) vector proportional to the external scalar and pseudoescalar 
fields. The constant B is related to the vacuum expectation value of quarks in the chiral 
limit and it can be expressed as B = m 2 /2m c where m is the pion mass at leading order 
in the chiral expansion. This notation follows that of Ref. 29) where the definition of the 
covariant derivative V M and the field strength tensor F^ v can be also found. Note that 
among all possible terms in £4 only those relevant for our calculation have been explicitly 
given. 

As a result of the procedure that leads to Eq. ( |25l) the explicit expressions for the LECs 
in terms of model parameters and form factors g p and f p can be obtained and eventually 
numerically calculated. However, a simpler and more direct way to obtain their numerical 
values is to extract them from the model predictions for the physical quantities they are 
related to. Of course, when proceeding in this way some care must be taken. In fact, since 
the Lagrangian in Eq. ( |25l) is valid up to fourth order in the chiral expansion, to extract 
the corresponding LECs from the physical quantities obtained in our quark scenarios we 
should treat them to the same order of approximation. Namely, we should calculate the 



relevant quantities in the chiral limit and its vicinity. As described in detailed in Ref. [22) . 
following this method the LECs t\, ...,£4 can be obtained from the predicted tt — n scattering 
parameters. Thus, in what follows we will concentrate on the remaining two parameters in 
Eq. f[2Tl) . namely £5 and £%. These are related to the pion charge radius (r 2 )^ and Fa (0) /m n 

by m 

B lff n F A (0) l /2 



4 V V2 



r 



'■) f 2 



^ = -l(r\L 2 - (28) 

The required chiral limit values of /„■, Fa (0) /m^ and (r 2 )^ for our three parameterizations 
of nonlocal quark models are given in Table III. The numerical results for £5 and £q are 
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given in Table IV, together with recent phenomenological values obtained in 3JJ and 



the values predicted within the Nambu-Jona Lasinio model [28[. For completeness we have 
also included the values of £i, ...,£4 from Ref. [22]. 

We observe that both the sign and the order of magnitude of the different LECs are in 
reasonably goo d ag reement with the results from %PT for \x ~ Im^. In fact, as already 



stressed in Ref. 



221 ] . for Set C the constants £2 and £4 are particularly very well reproduced, 



whereas the value of £3 is acceptable given the existing uncertainty. For £\ we obtain a 
rather large range of values when moving from one set of parameters to the other, but there 
is also a significative variation of this LEC with the choice of \x. Looking at the two new 
LECs, £5 and £q, we find that the (absolute) values obtained are relatively low, and the 
agreement with the phenomenological values from %PT can be achieved only for somewhat 
larger values of the renormalization point, fi > m p . This can be understood taking into 
account that our models do not include vector and axial-vector currents. Therefore, even if 
the quark couplings lead to some interaction in these channels, the models do not account for 
the effect of vector and axial-vector meson resonances. This is consistent with the analysis 
of Ref. 32| , which shows that £5 and £q are strongly dependent of the vector and axial- vector 
meson contributions. In fact, previous works on the pion electromagnetic radius including 
resonances indicate that vector meson contributions can give about 10—20% of the total 



result 



33, 



34]. 



V. CONCLUSIONS 

We have studied the vector and axial- vector form factors in 7r + — > e + v e r ) in the framework 
of nonlocal chiral quark models using three different sets of model parameters. As an 



extension of a previous work [23j where we have determined the values of Fy(0), -Fa(O) and 
Fy(0), here we have analyzed the vector form factor for values of q 2 up to 1 GeV 2 as well as 
the slope of FA{q 2 ) at threshold. While the latter is interesting in view of the comparison 
with future experimental results, the study of the behavior of Fy(q 2 ) presents relevant 
theoretical motivations, since it can be related to the form factor F 71 " 77 * (q 2 ) associated with 
the vertex 7177*. We have found that the model predictions for the g 2 -dependence of the 
vector form factor can be well described by a single-pole parameterization. This is indicated 
by the fact that there is a slight difference between the values of the parameters ay and Xy 
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for our three model parametrizations [the parameters ay and Ay are defined by Eqs. ( 123]) 
and 024p]. On the other hand, it is found that for all three parametrizations our results for 
ay and Ay are below experimental estimates. These discrepancies are not unexpected, since 
our model does not include vector- vector interactions. In fact, a simple estimate based on 
an extended NJL model [25| indicates that the vector contribution to 7r° — > 77* has the right 
order of magnitude and sign in order to account for the differences. This could be tested for 
the present nonlocal models that include wave function renormalization through the explicit 
inclusion of vector- vector interactions. Regarding the slope parameter A^, we find that 
our predictions are quite similar for the three parameter sets considered, and turn out to be 
about 3/5 smaller than the value obtained in a resonance effective model {4]. Finally, we also 
have extended the analysis of Ref . j^J , where the compatibility of our nonlocal models with 
effective meson theories such as %PT is studied. Indeed, we have been able to extract the 
model predictions for the low energy constants (LECs) £5 and taking into account 
the relation between these LECs and the values of -Fa(O) and the pion charge radius (r 2 ) n . 
Considering the results for £ 5 and t$ together with the values of £1, .., £4 previously obtained 
in Ref. 22], we observe that the predictions for both the sign and the order of magnitude of 
the LECs within the present non-local models are in overall reasonably good agreement with 
the results from ;^PT. For the particular case of £5 and we find that the (absolute) values 
obtained are somewhat small, a fact that as in the case of ay and Ay might be attributed 
to the lack of explicit vector and axial-vector interactions in our models. In any case, it is 
seen that nonlocal model predictions for £5 are closer to the result of ;\PT than the value 
obtained in the framework of the local NJL model. 
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APPENDIX I: Explicit expressions for ^(g 2 ) 



In this appendix we present the expressions for the contributions of the different diagrams 
in Fig. 1 to the axial form factor F^q 2 ). Let us start by the triangle diagram in Fig. la, 
which yields 

F A {q 2 ) /r~ f 



-N c V2G wqq J 



9ti 



z £+k z i-q z i 



{A b (A Cl + Bc 2 + Cc 3 ) 



(2tt) 4 D e+k D £ - q D e 

+B 5 (Ac 4 + Bc 5 + C c 6 ) + C 5 (A c 7 + B c 8 + C c 9 )} 



(29) 



where 



a = 1 

2 



B 



1 1 
+ — 

z e+k z e. 

1 



C 



2£ ■ k 
1 



k 



z l z l+k_ 

m e+k _ vfM 
z e+k z i 



A, 
C, 



1 1 
— + — 

z i z i- 



2q 



z l- 



1 

z t 



q 2 cr 2 
+ —a 



2 m q me m t - 



Zl-c, 



q 2 a 1 _ 



(£-,q) . (30) 



Here, and in what follows, we use the already introduced shorthand notation 
together with the functions a^(£,q) defined in Eq. f[2"Uj) . The coefficients q with i 
are given by 



±r/2, 



1,...,9 



c 1 = 2e 1 (me+k - me- q ) + 2 e 2 (me - me+k) - 2 e 3 (me + me- q ) + 4m e 
c 2 = e x [-2k ■ £ (2me + m £+k - me- q ) - 2d x + k ■ q (m e + m e+k )} 

+ 2 e 3 [d 3 - i ■ q (m e + me+k)] 
c 3 = ei [k ■ q - 2 k ■ £ + d A - d 6 ] - 2 e 3 [d 5 - 2£ 2 + I ■ q] 

c 4 = e 1 [-2di + 2£ ■ q (m i+k + 2m e - me- q ) - q 2 (m e + me+k) + k-q(m t - m £ _ 

+ 2e 2 [k ■ £ (mi^ q - me) + d 2 ] 
c 5 = ei [(Ak-£ + A£ 2 -U-q-k-q) d 1 - 2 k ■ £d 2 - 2£ ■ qd 3 

-2 k ■ ££ ■ q (Ame + m e + k - me- q ) + q 2 (2£ 2 (m t + me+k) + k ■ £ (3m £ + m e+k ))] 
06 = 6! [k-£[2d 5 -2£-q + q 2 ] + (2£ 2 - £ ■ q) d Q - d 4 (2£ 2 - U ■ q + q 2 ) - k ■ q d 5 ] 
c 7 = ei [d 6 - d 4 - 2 £ ■ q - k ■ q] + 2 e 2 [d 4 + 2£ 2 + k ■ £] 

c 8 = ei [2k-£(£-q-d 5 )- (2£ 2 + k-£) d 6 + (2£ 2 - 2£-q-k- q + k- £) d 4 ] 

c 9 = ei [di - k ■ £ (me- q - me) - £■ q{m e + m e+k ) - k ■ q me] , (31) 
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where 



d 1 = f (me + m e + k - m e - q ) + m e - q m e m e + k 
d 2 = I 2 (me- q + me +k - me) + m^, m t m e + k 
d 3 = £ 2 (me + me+k + me-q) ~ m t - q m t m e+k 



d 4 = m e me+k - d 2 

d 5 = m e me- q + I 2 (32) 

d 6 = me-q (m e + m e + k ) 



and 



ei 



6 q 2 (k ■ if 8k-££-q 2 i 2 

H — -; ; e 2 



2q 2 k-£ n i-q 2k- 
+ 2 - ; e 3 



(k ■ q) 3 (k • q) 2 k ■ q ' (k ■ q) 2 k ■ q 

The contributions of the other diagrams are the following: 



k ■ q 



(33) 



F A (q 2 ) 



m: 



F A (q 2 ) 



ml 



F A (q) 



F A (Q 2 ) 



m w e 



-^G -N 



dH ei 



(2tt) 4 £-k 



(2g £ -g e --g e +) 



z t + z e - m e + m r 



1 1 \ m f + — m 



z e+ z i 



lq f + I 2 



q <-q 

me m e+ m 



D e + D e - z f - z e + 



z i z e+ z e- 



m t+ m e - 
i-q 



= V2 G wg q N c 



0"i a 9 (4 9) • ^ + ^+ wi^- + a 2 a+ (£, g) t ■ £ q m e + - 1+ m e 
dH 



= V2 G nq q N c 



(2tt) 4 e-ke-q 
dH 



Zf-TUf- 

C k l k 

D - 



z„+m 



D 



9e+ + y a <!~(^) 



ei 



(2tt) 4 k-£q-£ 



Zl l me k z t m t 



X 



9i+ -9t--« ( a t ( £ ' 9) + a 9 ( £ ' ^)) 



-y/2G vqq N c 



dH \+ q \+ q 



(2-) 4 ^ ^ 



0"! 



k + q_ 



Vg (I, k, q) 
(k - qf 



-o 2 



m i- Z-Zk+q- m l+, 



^k + q ^ + 1 "^l+q " " k + 1. 



Vf k , g) 

(fc • g) 2 



1 

1 


1 


1 


'2k- 






t k + qJ 


1 




1 


1 


+ 2g- 






l q-k 



m e + (ei + e 2 - e 3 - 1) - m r (e x + e 2 + e 3 + 1) 



+ g 2 02 , g) 



fe + q 



m £ + (d - e 2 - e 3 + 1) 

fc + q 



- m r (ei + e 2 - e 3 - 1) > (34) 

k+q 1 
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where 

r] h {£,k,q) = {k.q) 2 [ dX [ dX' {e\ + e 2 X + e 3 X' + XX') h"+ (35) 

JO J-l \q + \'k 

for h = g or /. 
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TABLE I: Results for Fy(q 2 ) and its first and second derivative at q 2 = 0. All results should be 
multiplied by 1CP 2 . In column 5 we give the empirical values of m w + i 771 " 77 * (0) /\/2 and a' v . Note 
that in our model ay = a' v . In column 6 we give the prediction of the local Nambu-Jona-Lasinio 
(NJL) model. 





Set A 


Set B 


Set C 


Exp [2] 


Exp (vr 77*) [24] 


NJL 


Ref [4] 


F V (0) 


2.697 


2.693 


2.695 


2.58(17) 


2.80(8) 


2.441 


2.71 


ay 


1.91 


1.87 


1.98 




3.2(4) 






Ay 


1.651 


1.726 


2.011 


10(6) 




3.244 


4.1 


K 


0.020 


0.026 


0.046 











TABLE II: Results for Fa{o 2 ) and its first and second derivative at q 2 = 0. All results should be 
multiplied by 10 -2 . In column 5 we give the prediction of the local Nambu-Jona-Lasinio (NJL) 
model. 

Set A Set B Set C Exp [2] NJL From [4] 
F A (0) 1.319 1.614 1.825 1.17 ± 0.17 2.409 exp. input 
A A 1.22 1.17 1.26 - 1.97 

X' A 0.012 0.013 0.034 



TABLE III: Physical quantities used in the evaluation of q obtained in our different scenarios in 
the chiral limit. 

SA SB SC 
(r 2 )^ fm 2 0.335 0.325 0.316 
F A (0)/m n GeV- 1 0.097 0.118 0.133 
f n MeV 91.2 91.4 91.8 
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TABLE IV: Values of £{ obtained in our different scenarios. The values of l\ as a function of 



30, 



311 ] . The last two columns corresponds to the NJL predictions from 



u are obtained from Refs. 

n 

Ref. [28] for two different constituent quark mass: M = 220,264MeV 



LEC 


Non Local QM 




XPT 




NJL 




xlO 3 


SA SB 


SC 


fi = m p 


fJL = 2 


fi = m n 


M = 220 M 


= 264 


h 


-2.07 -1.39 


0.26 


-4.0 ±0.6 


-1.9 ±0.6 


-0.4 ±0.6 


-0.63 


-2.3 


4 


6.51 6.46 


6.41 


1.9 ±0.2 


6.2 ±0.2 


9.1 ±0.2 


6.3 


6.2 


4 


-1.1 -2.3 


-4.1 


1.5 ±4.0 


-1.8 ±4.0 


-4.0 ±4.0 


-8.5 


-3.5 


4 


15.0 17.2 


20.3 


6.2 ± 1.3 


19.1 ± 1.3 


27.9 ± 1.3 


22.7 


12.2 


4 


-4.39 -3.90 


-3.54 


- 5.22 ±0.06 


-6.29 ±0.06 


- 7.02 ± 0.06 


-2.88 


-2.60 


4 


-11.9 -11.6 


-11.4 


-13.1 ±0.4 


-15.3 ±0.4 


-16.8 ±0.4 


-11.5 


-10.4 


(2 4-4) 


3.15 3.82 


4.33 


2.7 ±0.4 


2.7 ±0.4 


2.7 ±0.4 


5.75 


5.21 
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(a) 




FIG. 1: Diagrammatic representation of the possible contributions to ir + — > e + v e ^ decay. Double- 
dashed lines, wavy lines and single-dashed lines represent the e v e pair, the outgoing photon and 
the decaying pion, respectively. While for the vector form factor only the contribution from the 
triangle diagram (a) is nonvanishing, in the case of the axial-vector form factor all five diagrams 
contribute. 
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q 2 [GeV 2 ] 



FIG. 2: F v (q 2 ) for sets A, B and C. The shadowed region corresponds to empirical data on F^ 11 , 
using Eqs. (J2TJ) and (32} 
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